Quantum Phase Transitions in the Shastry-Sutherland Model for SrCu 2 (B03)2 
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We investigate the quantum phase transitions in the frustrated antiferromagnetic Heisenberg model 
for SrCu2 (BC>3)2 by using the series expansion method. It is found that a novel spin-gap phase, 
which is adiabatically connected to the plaquette-singlet phase, exists between the dimer and the 
magnetically ordered phases known so far. When the ratio of the competing exchange couplings 
a(= J' /J) is varied, this spin-gap phase exhibits the first- (second-) order quantum phase transition 
to the dimer (the magnetically ordered) phase at the critical point a c i = 0.677(2) (a C 2 = 0.86(1)). 
Our results shed light on some controversial arguments about the nature of the quantum phase 
transitions in this model. 



PACS numbers: 75.10Jm, 75.40Cx 
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Two-dimensional (2D) antiferromagnetic quantum 
spin systems with the spin gap have been the subject 
of considerable interest. A typical compound found re- 
cently is SrCu2 (1303)2 PJ, in which the characteristic lat- 
tice structure of the Cu 2+ spins (see Fig. Q) stabilizes the 
singlet ground state. This system has been providing a 
variety of interesting phenomena such as the plateaus in 
the magnetization curve observed at 1/3, 1/4 and 1/8 of 
the full moment . The spin system may be described 
by the 2D Heisenberg model on the square lattice with 
some diagonal bonds which is referred to as the Shastry- 
Sutherland model ||, as pointed out by Miyahara and 
Ueda [Q. The key structure with the orthogonal dimers 
shown in Fig. [l] makes the system unique and partic- 
ularly interesting among 2D spin-gap compounds. In 
this frustrated system, there may occur non-trivial quan- 
tum phase transitions when the nearest-neighbor cou- 
pling J and the next-nearest-ncighbor coupling J' are 
varied. Albrecht and Mila || discussed the possibility 
of a helical phase between the dimer and the magneti- 
cally ordered phases by means of the Schwinger boson 
mean-field theory. Recent theoretical studies, however, 
have suggested that there may not be such a helical 
phase, but the first-order phase transition occurs from 
the dimer to the ordered phases ||J|. Furthermore, more 
recent study [Q claims that the phase transition should 
be of the second order with a non-trivial critical expo- 
nent v — 0.45(2). These controversial conclusions may 
come from the fact that the quantum phase transition 
in the Shastry-Sutherland model suffers from the strong 
frustration due to the competing exchange interactions 
J and J', and therefore a careful treatment should be 
necessary to figure out the correct nature of the phase 
transition. In particular, we have to keep in mind that 
such a strong frustration may possibly stabilize another 
spin-gap phase distinct from the dimer phase. 

In this paper, by calculating the ground state energy, 
the staggered susceptibility and the spin gap by means 
of the series expansion method, wc find that there should 



exist a novel spin-gap phase with the disordered ground 
state, which is stabilized by the strong frustration, be- 
tween the dimer and the magnetically ordered phases. 
The spin-gap phase found in this paper undergoes the 
first- (second-) order quantum phase transition to the 
dimer (the ordered) phase, when the exchange couplings 
J and J' are varied. The existence of the new phase 
can resolve controversial conclusions Q]||].[?J deduced for 
the quantum phase transitions in this frustrated model. 
We also point out that the material SrCu2 (1303)2 lies 
around the phase boundary between these two spin-gap 
phases, which may give a natural interpretation for the 
1/8-plateau formation in the magnetization curve. 

To investigate the frustrated spin system for the 
compound SrCu2(B03)2, we consider the 2D quan- 
tum Heisenberg model (Shastry-Sutherland model HU), 
which is described by the following Hamiltonian 



H — J ^ ' S,; • Sj I + J' ^ ' Sj • Sj, 



(1) 



where Si is the s = 1/2 spin operator at the i-th site 
and J (J 1 ) represents the nearest-neighbor (next-nearest- 
neighbor) antiferromagnetic exchange coupling. For later 



FIG. 1. Lattice structure of the Cu 2+ spins of 
SrCu2(B03)2- The nearest-neighbor bonds (J) are expressed 
by the solid lines and the next-nearest-neighbor bonds (J') 
by the dashed lines. 

convenience, we introduce the ratio a = J' /J. In Fig. 
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I], we have drawn the 2D Heisenberg model schemat- 
ically. We note that the system with only the next- 
nearest-neighbor coupling J' is equivalent to the Heisen- 
berg model on the square lattice which has a spontaneous 
staggered magnetization at T = [|||] . From this point 
of view, the nearest-neighbor coupling J is regarded as 
the coupling for a diagonal bond (see Fig. 0), which gives 
rise to the frustration together with J' ||]l| . 

In order to study the quantum phase transitions in 
this spin system, we employ the series expansion method 
developed by Singh, Gelfand and Huse p0| . We recall 
here that the quantum phase transitions in the Shastry- 
Sutherland model have been discussed by Weihong et al. 
H and Miiller-Hartmann et al. @| , by means of the dimer 
and the Ising expansions, from which the critical point 
between the dimer phase and the magnetically ordered 
phase has been estimated as a c — 0.691(6) and 0.697(2), 
respectively. As mentioned above, however, there is a 
controversy to be resolved about the nature of the phase 
transitions. Also, in order to determine the complete 
phase diagram, it is crucial to figure out whether there 
may exist another spin-gap phase besides the above two 
phases. We will address this problem in the following by 
using the series expansion method. 

To see our strategy clearly, we start with the 2D quan- 
tum spin model schematically shown in Fig. || HQ, 
which is topologically equivalent to the original model 
in Fig. [I]. In this figure, we have introduced an auxiliary 



FIG. 2. 2D spin system with the plaquette structure. The 
solid circle represents the s = 1/2 spin. The bold solid, the 
thin solid and the dashed lines represent the coupling con- 
stants J', XJ' and XJ 1 /a. When A = 1, this system is reduced 
to the Shastry-Sutherland model for SrCu2(BC>3)2. 

parameter A, which parameterizes the antiferromagnetic 
couplings labeled by the bold solid, the thin solid and the 
dashed lines, respectively, as J', A J' and XJ'/a(= XJ). 
Note that the system is reduced to the original Shastry- 
Sutherland model in the case of A = 1. An important 
point is that the introduction of A enables us to perform 
the cluster expansion starting from the isolated plaque- 
tte singlets (A = 0), which naturally describes the most 
likely spin-gap phase distinct from the dimer phase. 



To proceed the analysis based on the series expansion, 
we divide the original Hamiltonian eq. (|l|) into two parts 
as H = J' E S 4 • S 3 + \J2?ijSi ■ Sj], where T. l0 = 1 or 
a -1 for each bond on the square lattice (see Fig. ||). 
The first term is the unperturbed Hamiltonian which 
stabilizes the isolated plaquette singlets with the spin 
excitation gap. The perturbed Hamiltonian labeled by 
A connects these isolated plaquette singlets, by which a 
2D network develops. We expand the staggered suscep- 
tibility xaf, the spin-triplet excitation energy E(k) and 
the ground state energy E g as a power series in A. Here, 
to estimate the susceptibility, we introduce the Zeeman 
term H' = h \J2ieA &i ~ J2ieB &i] > where h is the stag- 
gered magnetic field and A(B) denotes one of the two 
sublattices. Note that an asymptotic analysis of the se- 
ries expansion is necessary to deduce the accurate phase 
boundary on which the susceptibility xaf diverges and 
the spin gap A = E(k = 0) vanishes. For this pur- 
pose, we make use of the Pade approximants Jll[] for both 
quantities obtained up to the finite order in A. Besides 
ordinary Dlog Pade approximants, we also employ biased 
Pade approximants JTT| , for which we assume that the 
phase transition in our 2D quantum spin models should 
belong to the universality class of the 3D classical Heisen- 
berg model H . Then the critical value of A c is determined 
by the formula xaf ~ (A c — A) 1 and A ~ (A c — A)" with 
the known exponents 7 = 1.4 and v = 0.71 Jl2| . 

We first calculate the staggered susceptibility xaf and 
the spin gap A by means of the plaquette expansion up 
to the fourth and the fifth order in A, respectively, for 
various values of a. Using the Dlog and the biased Pade 
approximants, we end up with the phase diagram shown 
in Fig. |^. In this figure, the solid (dashed) line repre- 
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FIG. 3. Phase diagram for the 2D spin system with the pla- 
quette structure in Fig. |^. The solid (dashed) line indicates 
the phase boundary obtained by biased Pade approximants 
for the spin gap (the staggered susceptibility). 

sents the phase boundary obtained by the biased Pade 
approximants for the spin gap (the staggered susceptibil- 
ity). When a — > 00 and A — 0, the system is reduced 
to an assembly of the isolated plaquettes with the spin 
gap. As A is increased, the correlation between these pla- 
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quettes grows up and the second-order quantum phase 
transition from the spin-gap phase to the magnetically 
ordered phase occurs at the critical point A c = 0.56 for 
a — > oo , which has already been studied by several groups 
p3| , P^Jl5) . On the other hand, decreasing a enhances 
the frustration, which in turn suppresses the antiferro- 
magnetic correlation, thus shifting the phase boundary 
upward for smaller a in the phase diagram. It is seen 
that two lines obtained from the distinct quantities are 
in good agreement with each other, which implies that 
the obtained phase boundary is rather accurate in spite 
of the lower-order pertubative calculation. By exploit- 
ing the phase boundary determined by means of biased 
Pade approximants for the spin gap, the critical value is 
given by a C 2 = 0.86(1) for A = 1. Recall that the sys- 
tem is reduced to the original model only for A = 1. We 
thus find that the Shastry-Sutherland model has the dis- 
ordered ground state in the region (0 < a < a C 2) on the 
A = 1 line. 

The above result does not necessarily imply that in 
the region < a < a C 2 the system always belongs to the 
disordered phase which is continuously connected to iso- 
lated plaquettes. In fact, it is known that the orthogonal 
dimer phase appears in the vicinity of a = |3|,|| . There- 
fore, it is necessary to clarify how these two spin-gap 
phases compete with each other by carefully comparing 
the ground state energy E s . To this end, performing the 
plaquette expansion up to the seventh order in A with a 
being fixed, we estimate the ground state energy E g for 
the Shastry-Sutherland model (A = 1) by means of the 
first order inhomogeneous differential method [ p~T[ j . The 
results are shown in Fig. |^. As mentioned above 
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FIG. 4. Ground state energy per site as a function of 
a = J'/J (A = 1, Shastry-Sutherland model). The flat line 
(E s / JN = —3/8) is the energy of the exact dimer state, while 
the solid line with dots (error bars are smaller than the line 
width) is obtained by the plaquette expansion. For compari- 
son, we also show the ground state energy obtained by Ising 
expansion Q as the dashed line. 

the system stabilizes the orthogonal dimer ground state 
for smaller a. It is found, however, that the first-order 
transition to the novel spin-gap phase introduced here 
occurs at the critical point a c i = 0.677(2). It is also seen 



from this figure that further increase of a induces the 
antiferromagnetic order, whose transition point is deter- 
mined by the crossing point of the ground-state energy 
obtained respectively by the Ising Q and plaquette ex- 
pansions. The result confirms the second-order phase 
transition deduced above, and the transition point es- 
timated from the figure is consistent with a C 2 = 0.86(1) 
obtained by the analysis of the susceptibility and the spin 
gap. Consequently, we end up with the phase diagram for 
the Shastry-Sutherland model as shown in Fig. ||. The 
present results shed light on the controversial arguments 
whether the quantum phase transition in this model is of 
the first or second order 0@^|. In those previous stud- 
ies, it was believed that the phase transition occurs only 
once between the dimer phase (I) and the ordered phase 
(III), giving rise to some confusions. Our phase diagram 
clearly resolves this problem by explicitly showing the ex- 
istence of the new spin-gap phase (II) which undergoes 
the first- (I^H) as well as the second-order transitions 
(II^III). 
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FIG. 5. Phase diagram for the Shastry-Sutherland model. 
The phase I represents the orthogonal dimer phase. The phase 
II newly obtained is adiabatically connected to the plaquette 
singlet phase. Ill is the magnetically ordered phase. 

To check the validity of the above phase diagram, we 
also show the results for the spin gap as a function of 
a = J'/J in Fig. @. In this figure, the results ob- 
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FIG. 6. The spin gap as a function of a = J'/J for the 
Shastry-Sutherland model. The solid line for a < a c i is the 
result obtained by Weihong et al. |j] The dots with error bars 
for a c i < a < a C 2 represent the spin gap at k = obtained 
by the plaquette expansion. 

tained by Weihong et al. || are shown for the orthogonal 
dimer phase (I: < a < a c i). In the new phase (II: 
a c i < a < a C 2), we determine the values of the spin gap 
at k = by means of the plaquette expansion up to the 
fifth order in A with the first order inhomogeneous differ- 
ential method. The results are shown as the dots with the 
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error bars. As seen in this figure, with the decrease of a 
from the second-order transition point a C 2, the spin-gap 
continuously grows up to stabilize the disordered ground 
state. As a is further decreased, the first-order phase 
transition occurs at a c \. 

In order to further confirm the present results, we have 
performed a different series expansion by choosing the 
isolated plaquettes with diagonal bonds as an initial con- 
figuration, which is different from the one shown in Fig. 
0. The calculation of the susceptibility up to the fourth 
order yields second-order transition with a C 2 = 0.87(3), 
being consistent with the above results. Furthermore, 
to confirm the first-order phase transition between the 
two spin-gap states, we have checked how the first-order 
phase transition point known for the ID orthogonal- 
dimer chain fl6| evolves with the increase of the inter- 
chain couplings. By performing the exact diagonalization 
studies for the 4x4 system, we have found that the first- 
order transition point for ID is continuously changed, 
and in the Shastry-Sutherland case, it coincides with the 
one found above within reasonable accuracy (a c i ~ 0.66), 
providing further support to our conclusion on the phase 
diagram. Although our results still seem to be partly 
contradicted to the staggered magnetization obtained by 
Weihong et al. ||, we believe that this could be resolved 
by further analysis of the results of the Ising expansion. 

Before concluding the paper, a brief comment is in or- 
der for the plateau-formation in the magnetization curve. 
Experimentally, the plateaus in the magnetization curve 
have been observed for the compound SrCu 2 (603)2 at 
1/3, 1/4 and 1/8 of the full moment In the theoreti- 
cal studies p7| , p8jl9 on the dimer phase, it has been clar- 
ified that the stripe order of the isolated dimer-triplets 
is important to understand the 1/3— and 1/4— plateaus. 
On the other hand, it is not so trivial why the 1 /8-plateau 
occurs in this compound, although a possible mechanism 
has been proposed |jl^,^8|. We think that the formation 
of the 1/8 plateau may reflect the fact that this com- 
pound is located around the first-order phase transition 
point between the two spin-gap phases and thereby pos- 
sesses the dual properties inherent in two distinct phases 
implicitly. We note here that the new spin-gap phase be- 
longs to the same phase as the Heisenberg model on the 
1 /5-depleted square lattice proposed for CaV40g. There- 
fore it is likely that the 1/8-plateau could occur in the 
same origin discussed by Momoi and Totsuka |Tj| for the 
plaquette system related to the 1 /5-depleted Heisenberg 
model. It is interesting to further clarify the mechanism 
of the 1 / 8-plateau by taking into account the above dual 
properties explicitly, which is now under consideration. 

In conclusion we have discussed the phase diagram 
for the Shastry-Sutherland model for the compound 
SrCu2(B0s)2 by means of the series expansion method. 
Our analysis has shown that there exists a novel spin-gap 
phase with the disordered ground state, which is adiabat- 
ically connected to the plaquette-singlet phase, between 



the dimer and the magnetically ordered phases known 
so far. When the exchange coupling ratio a = J'/ J is 
varied, the first-order phase transition occurs from the 
dimer state to the new spin-gap state, while the second- 
order phase transition occurs from this spin-gap state to 
the magnetically ordered state. This sheds light on the 
nature of the quantum phase transitions in this model, 
and resolves apparently controversial conclusions on this 
issue. 
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